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ABSTRACT 

In 2006, Prochter et al. reported a statistically significant enhancement of very strong Mg II 
absorption systems intervening the sightlines to gamma-ray bursts (GRBs) relative to the in- 
cidence of such absorption along quasar sightlines. This counterintuitive result, has inspired a 
diverse set of astrophysical explanations (e.g. dust, gravitational lensing) but none of these has 
obviously resolved the puzzle. Using the largest set of GRB afterglow spectra available, we reex- 
amine the purported enhancement. In an independent sample of GRB spectra with a survey path 
3 times larger than Prochter et al., we measure the incidence per unit redshift of > lA rest- frame 
equivalent width Mg II absorbers at z » 1 to be l{z)= 0.18 ± 0.06. This is fully consistent with 
current estimates for the incidence of such absorbers along quasar sightlines. Therefore, we do 
not confirm the original enhancement and suggest those results suffered from a statistical fluke. 
Signatures of the original result do remain in our full sample (£(z) shows an k 1.5 enhancement 
over ^(z)qso), but the statistical significance now lies at w 90% c.l. Restricting our analysis to the 
subset of high-resolution spectra of GRB afterglows (which overlaps substantially with Prochter 
et al.), we still reproduce a statistically significant enhancement of Mg II absorption. The reason 
for this excess, if real, is still unclear since there is no connection between the rapid afterglow 
follow-up process with echelle (or echellette) spectrographs and the detectability of strong Mg II 
doublets. Only a larger sample of such high-resolution data will shed some light on this matter. 

Subject headings: gamma-ray: burst - techniques: spectroscopic - quasars: absorption lines 
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1. Introduction 

In the last decade, the study of the inter- 
galactic (IGM) and circum-galactic medium (CGM) 
has received a great boost thanks to large spectro- 
scopic surveys of distant quasars, in particular the 
dataset provided by the Sloan Digital Sky Surveys 
( Yor k et al 1 l2000h . These objects randomly sam- 
ple thousands of lines of sight and, being bright 
background sources of light, probe gas and matter 
located in foreground objects. 

One of the most commonly surveyed set of tran- 
sitions in quasar spectra is the Mg II doublet at 
AA2796, 2803A. Its common detection stems from 
the large rest wavelength (which makes them eas- 
ily detectable by most optical spectrographs when 
the absorber is located at redshift z = 0.5 — 2.2), 
the relatively high abundance of Mg, and the 
strength of this resonance-line doublet. The Mg II 
systems are frequently classified in terms of the 
rest-frame equivalent width, W r , of the bluer 
component as "weak" (W2796 < 0-3 A), "strong" 
(^2796 > -3 A), and "very strong" (W 2 7Q fi > 1.0 
A, like in iRodriguez Hidalgo et al.l 120121) . For 
simplicity, throughout the paper we will refer to 
this last category as "strong" , since it is the only 
one pertinent to this work. Mg II doublet lines 
have been surveyed extensively from z w 0.1 — 2.5 
in the optical passband and now to z = 5 .2 



Qui der et alj2011t ISimcoe et al.l20lU IZhu fc Menard 



with near-IR spectroscopy (e.g. ISteidel fc Sargent 
19921 iNestor et all |2005i IProchter et al.l l2006at 



2012J). The results indicate that while the weak 
and strong absorbers incidence show small if any 
evolution with redshift, the very strong Mg II ab- 
sorbers present an increasing tren d up to z ~ 3 



before declining at higher redshi ft (jProchter et al 
2006at IMateiek fc Simcoel l2012[ ). This evolution 



rathe r closely tracks the cosmic star formation his - 
tory ijProchter et alJl2006at IZhu fc Menardll2012h . 



suggesting that some systems may b e causally con 



nected to on-going star format ion (jMenard et al 



201lt IMateiek fc Simcoel |2012|) , although, accu- 
rate analysis of the SDSS survey needs to be care- 
fully ta ken into account in o rder to avoid technical 
biases ( Lopez fc Chen 20121) . 

For several decades now, strong Mg II absorp- 
tion has been associated with gas in and around 
galaxies. Early work identified a small sample of 
L w L* galaxies at modest impact parameters 



(p ~ 10 — 50 kpc ) to quasars exhibiting strong 
Mg II absorp tion (Ber geron 11984 Lanzetta et al 



1987; Stcidcl 1993) , although no significant trend 
has been foun d for a population of Lumin ous Red 
Galaxies (e.g. iBowen fc Chelouchei 2011 , and ref- 
erences therein). 

These observations motivate the association of 
Mg II gas with the outer disk and/or CGM of 
these galaxies. Several QSOs line of sights pre- 
senting Mg 11 absorbers have been explored in 
order to probe the extent and the baryon con- 
tent around low-z ga laxies (jKacprzak et al. 2012; 
Chen fc Tinken 120081 . and reference therein), as 



well as a diagnostic of th e inner part of these galax 
ies' interstellar medium (|Bowen et al.lll995|). Also 



a stac k analysis was performed by IZibetti et al 



(2007) using light profiles (from associated galax- 
ies) of quasars exhibiting strong Mg II absorption 
in the SDSS. With their image-stacking technique 
they studied the cross-correlation between the 
Mg 11 gas and the galaxy light from 10 to 200 
kpc, finding that strong Mg 11 absorbers may be 
explained by models that include metal-enriched 
outflows from star-forming/bursting galaxies. 
Most recently, several attempts to trace the cover- 
ing fraction and nature of Mg II absorption by tar- 
geting known galaxies with coinc ident background 
quasars have been pe r formed ( Barton fc Cookd 
20091 : IChen et alJboid : IWerk et al.ll2012h . Their 



results indicate the mean covering fraction in- 
creases from - 70% for W r > 0.3 A to - 80% for 
> 0.1 A, confirming that extended Mg 11 absorb- 
ing haloes are a common feature around normal 
galaxies. Finally, it has been found that "strong" 
absorbers are often associated with nearby (within 
75 kp c) ~ 0.1 — 5L* galax ie s along the l ine of 
sight (iKacprzak et all l2008t INestor et ail l201lt 



Chen et al.|[201oh 



The survey and analysis of Mg II gas is no 
longer limited to quasar spectroscopy. For exam- 
ple, researchers have now used distant galaxies 
to probe foreground galaxies enablin g searches 
at v ery small impact parameter (| Rubin et al. 
20111 ) and statistical 'maps' of the absorption 



corre lated with the fore ground galaxy orienta- 



tion ( Bordoloi et al. 201 lh . Similarly, Gamma-ray 



Bursts (GRBs), with their extraordinarily bright 
optical afterglows provide not only direct infor- 
mation on their host galaxies, but also trace mat- 
ter intercepting their lines of sight ([Metzger et al 
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19971) . 

The advantage of using GRBs as background 
sources is twofold: first, the y can be observed 
up to very high redshifts (Kawai et al. 2006; 



along QSO lines of sight (iMenard et al.l [200p 



IPorciani et alll2007at iBudzvnski fc Hewettll201ll ): 



3) a geometric effect difference due to the sizes 
of the emitting regions between GRBs and QSOs 



Tanvir et al.l200alSalvaterra et al.l200alCuccniara et al.l (IFrank et al.l200ft IPorciani et alJ2007aHLawther et al 



2011bl ). which allows one to explore a larger red- 
shift path length, and second, their discovery is 
largely unbiased with respect to intrinsic proper- 
ties of their hosts (extinction, luminosity or mass) . 
When a GRB fades away, they leave the line of 
sights clear for future deep observ ations in order to 



search for the Mg II counterparts jV rccswiik et al 



2003; 
Chenl 



Jakobsson et al.l 12004 ISchulze et al.l l2012t 
20121 ). One of the first attempt to iden- 



tify the nature of three absorbers alon g GRB 
060418 was performed by IPollack et ail (|2009h . 
which identified the absorbers to be L ~ 0.1 — 1L* 
galaxies at very small impact parameter from the 
GRB location (p < 10 h^ 1 kpc). Deep imaging 
of severa l other fields have confirmed these early 
findings (|Chen et al. 2009). On the other hand, 



the number of GRBs discovered and spectroscop- 
ically observed is several orders of magnitude less 
then the number of quasars available in large opti- 
cal surveys (e.g. SDSS DR8). This difference has 
been reduced, however, with the success of the 
Swift satellite providing the d iscovery and follow- 
up of several hundred GRBs ( Gehrels et al.l 2004L 
l2009h . 

Shortly after the launch of Swift in November 
2004, a survey of Mg II absorption in GRB af- 
terglow spectroscopy was performed for an early 
sample of Swift bursts and a hete rogeneous sam- 
ple of pre-existing GRB spectra ( Prochter et al.l 
2006b, P06 hereafter). The authors revealed an 
extremely puzzling result: the incidence of strong 
(W2796 > lA) intervening Mg II absorbers was 
about 4 times higher along GRB sightlines than 
quasar sightlines. Despite the small sample size, 
the statistical significance of their dataset was 
high: the null hypothesis that GRBs and quasar 
spectra would show identical incidences of strong, 
foreground Mg II absorption was ruled out at 
> 99.99% confidence. 

The authors proposed several hypotheses that 
might explain the difference, which have since 
been studied in greater detail: 1) a possible in- 
trinsic origin of these ab sorbers associated near 



the GRBs themselv es ( Cucchiara et al. 1200a" 



2012 ) ; 4) a gravitational lensing e ffect (IVergani et al.l 



Porciani k, Madau 2001 : iRapoport et al 



2009; 

120111 ). Subsequent work has ruled out several 
of these possibilities and none appears to be suffi- 
cient on its own to explain the observations. After 
seven years of the Swift mission and more than two 
hundred GRBs with spectroscopic confirmations, 
this mystery remains. 

It is important to emphasize that the original 
P06 work, and even the studies that have followed, 
have relied on a small sample of GRB afterglow 
spectra. Even the largest analysis to date analyzed 
only 26 line of sights (finding 22 ab sorbers), for a 
total redshift path of Az — 31.55 ( Vergani et al.l 
2009). Furthermore, no study has analyzed a com- 
pletely independent set of GRB sightlines from the 
P06 analysis. 

In this paper, we use data obtained primarily 
during the Swift era by several facilities, to obtain 
the most complete sample of GRB afterglow spec- 
tra and the largest redshift pathlength available to 
date. From this parent sample, we are able to con- 
struct sub-samples which are entirely independent 
from the original work of P06. Similarly, we can 
study possible instrumental biases (e.g. spectral 
resolution) which may affect the final results. 

The paper is structured as follow: in Sj2]we de- 
scribe our dataset and the data analysis procedure, 
while in £j3] we present our procedure for defin- 
ing the redshift path density per GRB sightline. 
21 describes the search methodology to identify 
possible Mg 11 systems along every line of sight, 
with distinction between different datasets (e.g. 
high-resolution vs. low-resolution, strong vs weak 
Mg 11 equivalent width). Finally, in Section [J5] 
and Sjniwc present our findings, including interest- 
ing sub-samples results, and we summarize them 
in light of possible steps forward into understand- 
ing this puzzling phenomenon. All the quoted er- 
rors, unless otherwise stated, are considered at 1-a 
confidence level. 



Bergeron et al.l 120111 ): 2) a significant dust bias 
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2. Data Selection 

The acquisition of an optical spectrum from a 
given GRB afterglow is a complex and unrepeat- 
able process. During the Swift era, the timelapse 
between the discovery of the gamma-ray emission 
(by the Swift /BAT instrument) and the afterglow 
localization (by XRT and/or UVOT on-board the 
spacecraft) is generally less than a few minutes, 
with some exceptions due to observability con- 
straints which delay the satellite to slew towards 
the BAT position (e.g., due to the small angular 
separation between the GRB and the Moon or the 
Sun). 

The on-board localization has an accuracy be- 
tween several arcminutes (BAT only) to subarc- 
seconds (UVOT) . The immediate transmission to 
the ground of all the Swift-ac quired data via 



of the GRB redshift. For this purpose, the identi- 
fication of fine-structure lines represents a secure 
determination of the GRBs host galaxies and the 



the G amma- Ray Burst Network (jBarthelmv et al 



19951 ) allows rapid (seconds to hours) follow-up 
with ground-based optical/IR telescopes. 

In most of the cases presented in this paper, 
rapid follow-up spectroscopic observations were 
triggered as soon as an X-ray counterpart posi- 
tion was delivered (XRT identified more than 98% 
of the BAT GRBs). Before the actual trigger is 
sent to the opportune telescope, we check the cur- 
rent telescope/camera set-up available, the visibil- 
ity window at the telescope site, and weather con- 
ditions. A finding chart is usually provided to the 
telescope operator, using archival images (usually 
SDSS or USNO catalogues), which also helps to 
identify the afterglow. Once the trigger is sent to 
the telescope, usually via Target of Opportunity 
(ToO) programs, a long acquisition image of the 
field is acquired. In the eventuality that a "new" 
source has been found inside the XRT error circle 
(by comparison with the finding chart), a spectro- 
scopic sequence is executed. 

In other cases, especially for high-z bursts, 
robotic, real-time follow-up by different facili- 
ties have provided similarly accurate identification 
within the first hour, using redder filters than the 
ones available on Swift /UVO T (e.g. the GROND 
and RAPTOR instr uments, iGreiner et ail 120081 : 
Vestrand et al. l2002h . Thanks to the prompt re- 



sponses, different groups have been able to obtain 
spectroscopic observations of the optical afterglow 
when it was bright enough to detect absorption 
lines, which generally yields a definitive estimate 



GRB circumstellar environment (jProchaska et al 
20061) . Other secure identification is the pres- 



ence of a Damped Lyman-a system, which also 
has been signature of a typical high-z GRB host 
galaxy environment. Whenever these features are 
not present, we assume that the higher redshift 
system of absorption features from several ion- 
ized transitions corresponds to the GRB redshift, 
but obviously does not guarantee that these fea- 
tures do not rise, instead, in foreground objects 
along t he GRB line of sigh t (see for example GRB 
071003. IPerlev et al.ll2008h . 

For GRBs that occurred before the launch of 
the Swift satellite a similar procedure was fol- 
lowed: the gamma-ray identification was made by 
high-energy facilities in orbit (BATSE, HETE-2 
spacecrafts), while X-ray observations were due to 
slower response missions (like Beppo-SAX) , which 
required several hours to repoint. This inevitable 
delay propagated often to a late follow-up obser- 
vation, which, in few cases, led to a spectroscopic 
sequence to executed. 

For our parent sample, we set out to obtain 
the optical afterglow spectra for all GRBs with re- 
ported redshifts, restricted as follows. Because of 
our interest in detecting Mg n lines, and because 
most spectrographs have wavelength coverage be- 
ginning at ~ 4000 A (or poor UV sensitivity), we 
require the Mg n doublet rest-frame wavelength to 
be redshifted beyond this limit. This leads us to 
include all the publicly available spectra obtained 
from GRBs with redshift higher than zqrb = 0.5. 

The spectra analyzed in this paper were ob- 
tained with facilities across the world, including 
the Gemini Observatory, Keck Observatory, and 
Very Large Telescope. Many of these data were 
obtained by our respective research teams, al- 
though several tens were taken from public data 
archives or were kindly contributed by members 
of the community. Table [T] lists all of the sight- 
lines with reported GRB redshifts where we were 
able to retrieve a spectrum. The last column lists 
the literature references for GRB afterglow spectra 
that have been previously published. 

A small sample of seven GRB spectra, as men- 
tioned in was obtained during the pre- Swift 
era: these GRBs were discovered by non-GRB 
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dedicated missions, like the Interplanetary Net- 
work (GRB 000926), Beppo-SAX (GRB 010222) 
and HETE-2 (GRB 020813,GRB 021004,GRB 
030226,GRB 030323,GRB 030429), and followed 
from the ground several hours (if not days) af- 
ter the events were discovered. Neverthless, these 
data have sufficiently high quality to be included 
in our work. 

Our large dataset consists of a total of 118 GRB 
afterglows observed by different facilities and in- 
struments, including the Gemini telescopes with 
the GMOS instruments (46 spectra), the Very 
Large Telescope equipped with FORS1, FORS2, 
UVES, and the X-Shooter spectrographs (55), the 
North Optical Telescope (NOT) with the AL- 
FOSC camera (6), the Keck telescopes with the 
HIRES,LRIS and ESI instruments (13), and the 
Magellan Clay telescope with the MagE spectro- 
graph (1). Finally, we also include two spec- 
tra obtained by the KAST spectrograph mounted 
on the 3-m Shane Telescope at Lick Observatory 
and several from the Magellan and the Telesco- 
pio Nazionale Galileo (TNG). The spectral reso- 
lution of these data ranges from 450 km/sec (or 
- 13A , NOT/ALFOSC) to 7 km/sec (~ 0.13A, 
HIRES/UVES). This large variety of data give us 
the opportunity to test different subsamples drawn 
from the overall 118 GRBs. All the data presented 
are part of a public repository of GRB spectre^]. 

The VLT sample has been obtained with a spe- 
cific set of observational criteria, making the best 
use of the different instruments available and the 
technical improvement for ToO observations, like 
the Rapid Response Mode (RRM), which allows 
the observer to execute the required observation 
remotely and without a major intervention by the 



telescope operators (see Fvnbo et al. 2009() 



Most of the FORS 1/2 data are part o f the 



catalog presented in iFynbo et al. (|2009h and 



de Ugarte Postigo et al.1 ~( 2012a[ ). while most of 
the high-resolution ones (HIRES, ESI, UVES) 
were al ready published in sin gle-GRB papers or as 



part of IVergani et al 



(|2009f ). A large fraction of 
the Gemini spectra are presented here for the first 
time, and are the result of our group' s follow-up 



Cucchiara 



effort s over the last 7 years (see also 
2010l ). We encourage the readers to refer to the 



reference in the last column in Table[T]for the data 



reduction procedures and the original published 
papers. In the following sections we will briefly 
review the reduction procedure for the Gemini 
and the X-Shooter data. 

2.1. Gemini sample 

These datasets are part of several follow-up pro- 
grams for which Target of Opportunity time was 
awarded between 2005 and 2011. All the data 
included were obtained with the Gemini Multi - 
object spectrographs (GMOS: iHook et alJliool . 
The typical observation sequence consists of two 
spectra in two dithered positions along the slit 
(usually 1" wide) in order to facilitate sky-line sub- 
traction. Immediately before or after the science 
frames, a ThAr lamp is observed and a flat field is 
obtained in order to allow data reduction "on-the 
fly". 

We used the GEMINI /GMOS data analysis pack- 
ages under the iraeH environment in order to per- 
form the basic reduction, flat fielding and wave- 
length calibration. Cosmic rays were identified 
and replaced by a median of the surrounding pixels 
which were not flagged as bad pixels. For this pur- 
pose we used the lacos_spec tool (jvan Dokkum 
2001). Finally, one frame was subtracted from the 
other to remove the strongest skylines. This pro- 
cedure provides good results at A < 8000A, but 
leaves significant residuals in the reddest portion 
of the spectra where the GMOS spectrographs suf- 
fer substantial CCD-fringing. Therefore, the ex- 
tracted error arrays associated with the Gemini- 
GMOS data reflect these higher-noise patterns at 
longer wavelengths. 

One dimensional spectra were then extracted 
using the IRAF APALL tool and coadded weight- 
ing each spectrum by the inverse of its variance 
spectrum in order to increase the S/N of the fi- 
nal result. The APALL package also produces a 
1-d array with the poissonian statistical error and 
the 1-d sky background (estimated in regions se- 
lected far from the object trace, so to avoid any 
spurious contamination). These last two arrays 
have been summed in quadrature to obtain the 
final error array per pixel. In some cases we as- 



: http : //grbspecdb . ucolick . org/ 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association for Re- 
search in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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sess the quality of the extracted error array with 
the estimated RMS of the data-array and modified 
the latter in order to fully account for the poisso- 
nian fluctuations in the actual data. Finally, us- 
ing the splot routine we estimated signal-to-noise 
over the whole wavelength range (also reported in 
Tabled]). 

2.2. X- Shooter data 

Data for GRB 090926A and GRB 100418A, 
were obtained via the ESO Archive! and re- 
duced with version 1 .3.7 of the X-Shooter pipeline 
dGoldoni et al.ll2006h in physical mode. The spec- 
tra in the UVB and VIS arms were used for the 
redshift pathlength estimate as well as the Mg II 
search. We do not use NIR arm due to the high 
level of contamination from skylines in the in- 
frared. Furthermore, the infra red sample of QSO 
spectr a, largely obtained by iMateiek fc Simcoe 
(|2012l ). are still small compared to the large com- 
pilation from SDSS. 

2.3. Subsamples 

The set of spectra listed in Table Q] comprises 
our full sample for analysis which we refer to as 
Sample F. This sample maximizes the survey path 
for Mg II absorption along GRB sightlincs. From 
this parent sample, we consider several subsamples 
for the same analysis. Most important is the in- 
dependent subsample (Sample I) which ignores all 
of the data analyzed in the original paper of P06. 
We focus first and foremost on this subsample to 
perform a complementary study. In addition, we 
consider two other subsamples which cut the data 
according to spectral resolution: we combined all 
the high-resolution spectra, obtained with echelle 
or echellette spectrographs (ESI, HIRES, MagE 
and UVES) in Sample H and, all other data in 
Sample L. These are summarized in Tabled 

3. Survey Path 

The starting point of a survey for interven- 
ing absorption-line systems is to estimate the red- 
shift path density g{z). This function expresses, 
as a function of redshift, the number of unique 



sightlines for which an absorption line could be 
detected in the survey to a limiting equivalent 
width. In practice, one determines for each spec- 
trum those regions that have sufficient S/N and 
are free of strong blending by terrestial or intrin- 
sic gas. These specific windows define redshift in- 



tervals j, 



-Xl ^2\3i 



for the ith quasar (or GRB) 



3 Based on observations made with ESO Telescopes at the 
La Silla Paranal Observatory under programs ID 60. A- 
9427(A) and ID 085.A-0009(B) 



where g(z) = 1 within each window and zero oth- 
erwise. By integrating g(z) across the full spec- 
trum, one recovers the redshift path Azi covered 
by the source. 

To properly determine g{z) for each sightline, 
several issues must be considered to minimize sys- 
tematic effects that could bias the search. First, 
we exclude from the search the wavelengths (or 
the redshift ranges) that fall in the atmospheric 
telluric bands, which heavily absorb the afterglow 
flux, rendering it very difficult to identify any fea- 
tures (intrinsic to the GRB, or QSO, host or inter- 
vening). Since some of our spectra extend towards 
the near-IR regime, we also consider atmospheric 
absorption at these wavelengths. A complete list 
of the excluded regions is presented in Table [3] 

Second, GRB afterglow spectra exhibit strong 
absorption lines belonging to ionic species located 
in the progenitor environment and up to tens kpc 
along the line of sight. The number of detected 
host features varies depending on the brightness 
of the afterglow, the properties of the host galaxy, 
the signal-to-noise and t he resolution of th e spec- 
trograph (see Figure [T]). IChristensen et aL ( 2011 ) 
created a high-S/N composite spectrum using 66 
afterglow spectra obtained with low and mid- 
resolution spectrographs. Strong absorption lines 
were identified as well as weak ones previously un- 
detected in the individual spectra. Since most of 
these lines are common in GRB host galaxies we 
compile a sublist of these absorption features to 
be excluded in our redshift pathlength calculation. 
We included also some of the most common fine- 
structure transitions. To be more conservative, 
a region equivalent to one half-resolution element 
both blueward and redward of the observed cen- 
tral wavelength of the considered transition (as set 
by the redshift of the host galaxy) has had g(z) set 
to zero. In the cases of high-resolution spectra, the 
minimum size of the masked region is 200 km/s. 
A complete list of these features is also presented 
in Table El 

Finally, for some of the spectrographs (e.g. 
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UVES, GMOS), the spectral coverage is non- 
contiguous due to gaps between detector chips 
and/or the use of multiple cameras. Regions with- 
out data were simply masked in the g{z) evalu- 
ation. In several cases, regions beyond - 8000A 
were heavily effected by fringing, even after cor- 
recting for it in the data processing. We opted for 
a visual inspection of the data and decided on a 
case-by-case which regions needed to be excluded 
for the search. 

As an example, in Figure[2]we present the Gem- 
ini/GMOS spectrum of GRB 060210, where all 
of the masked regions are indicated. It is clear 
from Figure [2] that the maximum redshift, z max , 
allowed for our intervening Mg 11 search is dic- 
tated by the host galaxy redshift, in particular, 
we begun our search starting 1500 km s _1 blue- 
ward the corresponding Mg 11 feature (or z max — 
zgrb — 0.015). Also, the minimum, z m i n , is in- 
dicated either by the bluest wavelength covered 
by the spectrograph or, as in the case of GRB 
060210, by the presence of the Lya feature (at 
rest-frame Affy* = 1215.67 A) such that z min = 



(if-l) + 0.05, where Xf y s a = Xlf a (1 + z GRB ) 
(or equivalently ~ 5000 km s _1 redward the Lya 
feature). We do not extend the search for inter- 
vening Mg II into the Lya forest and we avoid 
the (typically) very strong damped Lya absorp- 
tion profile of the GRB host galaxy. 

Once these regions have been excluded we de- 
termined the 5cr equivalent width limit per pixel, 
using the variance spectrum associated to each ob- 
ject (shown as red in Figure [2]), considering a sim- 
ple gaussian profile of FWHM — 5/2.35, where 
5 is the resolution element. At each unmasked 
pixel in the spectrum, we query whether the 5a 
equivalent width limit exceeds a given rest-frame 
survey limit for Mg II 2796 (e.g. lA). If the limit 
is satisfied, we query whether the corresponding 
Mg II 2803 line lies in an unmasked region. If both 
of these criteria are satisfied, g(z) = 1 for the red- 
shift interval covered by that pixel otherwise we 
set g(z) = 0. This generally leads to a series of 
discontinuous redshift intervals for the Mg II sur- 
vey, as listed in Table 0] 

Figure [3] presents the total redshift path den- 
sity for Sample F and Sample I, which represents 
the number of GRB sight-lines available for our 
Mg 11 search as a function of redshift. These 



are shown for a limiting rest-frame equivalent 
width of lA at 5cr confidence. It is immediately 
clear from this figure that we accurately excluded 
from our analysis the telluric lines regions (i.e. at 
ZMgii ~ 1-5, 1.7, 1.9, and 2.4). Also, the analysis is 
mainly performed where we have the majority of 
the searchable path, in the 0.4 < z < 2.2 interval 
range, due to the larger statistical sample from the 
GRB and the QSO samples. 

The total Az of the survey crudely expresses 
its statistical power. This may be calculated by 
simply summing the Azi values for each source. 
For the full sample (F), a redshift pathlength of 
Az = 55.5 for the lA equivalent width limit. For 
the independent sample (I), we find Az = 44.9. 
The latter represents a ^4 times larger survey 
sample than P06. 

4. Identifying and measuring intervening 
Mg 11 absorbers 

We described in the previous section the con- 
struction of the redshift path density which de- 
fines, for each spectrum, the regions where an in- 
tervening Mg 11 doublet may be detected at 5a 
significance. Independent of this calculation, we 
have searched each sightline for the presence of 
Mg 11 absorbers. 

Using the algorithm for optimal extraction 
(|Hornelll986[ ). we constructed an equivalent width 
spectrum by convolving the normalized data with 
a Gaussian profile with width set by the resolution 
of the spectrograph and weighting the flux at each 
wavelength by the associated variance. From this 
array, we determined every feature satisfying a 5a 
detection threshold via an automatic procedure 
similar to the C IV doublet search performed in 
Cooksev et al. ( 2010f ). We considered each line as 
a possible Mg 11 absorber, and confirmed this as- 
sociation through the presence of a proper Mg 11 
doublet (both in velocity separation and relative 
W r ). Finally, we inspected every candidate iden- 
tified by this procedure visually, confirming the 
presence of a genuine doublet with the additional 
identification of other common features (e.g. Mg 1 
Fe 11). We also accurately measured the equiva- 
lent widths of the doublet components via line 
profile fitting (see Table [5]). As sanity check, each 
GRB sightline was also manually inspected by 
the lead authors in search of Mg 11 doublets that 
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might be missed by the automatic screening pro- 
cess. We found two doublets in addition to the 
candidates automatically identified which may be 
Mg II features (see tH.ll for our completeness anal- 
ysis). Also, we found some doublets which were 
misidentified as Mg II doublets: in reality these 
features were host galaxy fine-structure transitions 
or other metal lines belonging to other interven- 
ing systems (e.g. GRB 061121). We estimated 
that our total redshift pathlength would be de- 
creased of a factor of < 6% if we would have 
masked also these features, therefore we prefer 
not to exclude these spectral regions to preserve 
a maximum searchable path. Again, our visual 
inspection prevent these features to be accounted 
in our Mg II search. 

Table [3] lists the Mg n systems that have been 
discovered, and Figures llOaHlOml show the line- 
profiles of all the strong Mg n systems in combi- 
nation, when available, with other metal features. 
The W r values for the Mg n doublet were esti- 
mated by fitting the line-profiles with Gaussian 
profiles or in the case of line-black saturated tran- 
sitions (e.g. high resolution data) by pixel sum- 
mation. We estimated the uncertainty in our W r 
values by summing the pixel-by-pixel variance in 
quadrature. 

4.1. Completeness estimate 

It is important to note that the sample un- 
der consideration has been obtained by a large 
variety of facilities and, also, GRBs have been 
observed at different epochs (meaning at differ- 
ent afterglow brightness) as well as with different 
atmospheric conditions. Therefore, it is worth- 
while evaluating our completeness in finding very 
strong Mg II absorbers at the considered 5er con- 
fidence level. To assess our completeness we in- 
serted mock Mg II features into our spectra (tak- 
ing into account the S/N and the resolution of the 
original spectra) and then we re-process these new 
datasets via our automatic procedure. The in- 
jected features, in a number which is drawn by 
a poisson distribution centered on the expected 
number of absorbers (£(z)xAz) for the GRB sam- 
ple, have random equivalent widths between 0.05 
to 5A and a maximum number of seven sub- 
components, each with a range of doppler param- 
eters b — 5 — 20 km/s. These features were in- 
serted between z min and z max as defined in SecJ3] 



per each GRB. We repeated this process 50 times 
per sightline for a total of 5250 iterations. We 
compared the number of injected strong features 
(W r > lA) that should be automatically identified 
because they were located in regions of the spec- 
tra were g{z) — 1, (accordingly to Sec. [3]) with 
the actual recovered list: we conclude that « 98% 
of the systems were correctly identified and de- 
tected as genuine strong Mg n doublets. Figure 
@] shows the result of our completeness test: on 
the top panel we show the total number of ab- 
sorbers correctly identified (black histogram) and 
not (in red) depending on the instrument resolu- 
tion. It is clear that the lowest-resolution spectra 
(e.g. the V300 grating with the FORS1 spectro- 
graph) have a low completeness level using our 
automated search but that the other spectra give 
excellent results. For the missed absorbers in the 
lowest-resolution data, we find that these doublets 
are usually self-blended (resembling a single broad 
line) or they are blended with the profile wings of 
other lines (intervening systems metal lines or host 
galaxy features) preventing the automatic identi- 
fication of both doublet components. 

We examined again our original sample and we 
confidently retrieved only two such cases: GRB 
090812 and a possible absorber at z = 1.055 
and GRB 070110 with a possible doublet at z = 
1.5875. Nevertheless, including such features, 
which were not automatically recovered, does not 
effect our conclusions. 

In the bottom panel of Figure [4] we present our 
cumulative completeness level with increasing res- 
olution. Again, as noticed previously, we reach 
~ 97% level around the resolution of the Gemini- 
GMOS instrument (R400 grating, R ~ 1200), 
whose spectra provide the best combination of 
Signal-to-Noise ratio and resolving power to prop- 
erly identify the population of strong Mg II ab- 
sorbers characterizing the GRB intervening sys- 
tem population. 

5. Results 

5.1. Incidence £(z) 

Combining the results from the previous two 
sections, we may estimate the incidence of Mg n 
absorption per unit redshift £(z) (also referred to 
as dN/dz or dn/dz). The standard estimator for 
£(z) to a limiting W r is the observed ratio of the 



S 



number of absorbers discovered, N, having W > 
W r in a given redshift interval [zi, Z2] to the total 
redshift pathlength searched, Az, in that redshift 
interval 



t{z) 



with 



Az 



N 
Az 



g{z) dz 



(1) 



(2) 



Figure [5] presents our £(z) estimates for W r > lA 
Mg 11 absorbers, for Sample I and Sample F re- 
stricted between z = 0.4 and 2. The error esti- 
mates assume Poisson statistics for N and corre- 
spond to 68% confidence. The values for the GRB 
sightlines have roughly constant value with red- 
shift at £(z)grb ~ 0.18 and £(z)grb ~ 0.36 for 
Sample I and Sample F respectively. 

For comparison, we display a fit to the mea- 
sured £(z)qso values for W r > lA Mg 11 absorbers 
discovered along the thous ands of quasar sight - 
lines drawn from the SDSS ( Zhu fc Menardl l2012). 
These spectra have been chosen to have S/N > 15, 
so to assure a high confident statistical sample of 
Mg 11. 

This quasar sample was searched for Mg II ab- 
sorbers at wavelengths redward of the strong C iv 
quasar line (A rest = 1550A) and blueward of the 
reliable response of the Sloan fibers up to the 
quasar's Mg II emission line. 

From Figure [3] it is evident that while the Sam- 
ple I follows the expected distribution derived 
from the QSO analysis, Sample F still presents 
a modest excess of absorbers. In the case of Sam- 
ple F, for instance, Az = 55.5 and the number 
of absorbers identified is N b s = 20 (N exp = 13). 
Overall the €(z)grb,f = 0.36 ±0.09, a factor - 1.5 
greater than the expected quasar density of ab- 
sorbers (£(z)qso,f = 0.24). Considering the inde- 
pendent sample, which as mentioned in §2.31 ex- 
cludes all the lines of sight in P06, we obtain 
^(z)grb,i = 0.18 ± 0.06. Following the same anal- 
ysis, similar results are evident using the high- 
resolution and the low-resolution samples (Sam- 
ple H and Sample L): in these cases we identify an 
overabundance of strong Mg 11 absorbers in the 
high-resolution sample, leading to a £(z)h = 0.64, 



a factor 2.6 larger then the expected (^(z)qso,h = 
0.25). We summarize our analysis in Table [2] 

Figure [6] shows the cumulative distribution of 
Mg 11 absorbers detected from GRB Sample F and 
Sample I, together with the quasar estimates. We 
may compare these results against the predicted 
cumulative distribution functions for a QSO sur- 
vey with identical search path to the GRB analysis 
by simply convolving the GRB g(z) with £(z)qso- 



N 



QSO 
cumul V 



>z')= / e{z) QSO g{z)dz . (3) 
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It is evident that the full Sample F exhibits a mod- 
est excess of ~ 30%, but that the independent 
Sample I shows no excess. The new results for 
Sample I do not confirm earlier works which re- 
ported an excess of strong Mg II absorption along 
GRB sightlines. 

5.2. Monte Carlo Analysis 

To assess the significance of these results, in 
particular the observed excess for Sample F, we 
perform a Monte Carlo analysis as follows. First, 
we selected a s et of 12700 SDSS quasars from 
Zhu fc Menardl (|2012h that have a continuous 
g(z) — 1 redshift path density from z m ; n to z max , 
where z m \ n is the greater of 0.4 and (1 + zqso) x 
AgfvAMgii and z max = min[z QSO - 0.04,2.2]. 
This is the brighter subset of quasars in the SDSS 
with correspondingly higher S/N spectra. Re- 
stricting our Monte Carlo analysis to this QSO 
sample facilitates the generation of random sam- 
ples with a survey path identical to the GRB 
analysis. 

For GRBs with z < 1.5, an SDSS QSO matched 
in redshift will cover the survey path of the GRB 
analysis. For a given GRB, we selected all quasars 
close in redshift space to zqbb (usually in the 
range z GRB < zq SO < z GRB + 0.04 there were 
always at least 50 such quasars). In each Monte 
Carlo realization we randomly picked one and 
by construction adopted the g(z) from the refer- 
ence GRB spectrum. We then ide ntified the total 
numb er of absorbers discovered bv lZhu fc Menardl 
(|2012l ) along the lines of sight of these quasars and 
recorded those that satisfy the W r > lA limit and 
have g(z) = 1. 
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For zqrr > 1-5. the Mg II survey performed 
bv lZhu k Menardl f)2012h using the SDSS quasars 
does not extend as low in redshift as our GRB 
analysis because those authors truncated the 
search bluer then the C IV emission peak. As 
a result, we considered two approaches to han- 
dling this difference. The cleanest approach is to 
artificially truncate the GRB analysis at the same 
starting redshift as the quasars, i.e., 



z min,GRB 



(zqso + l)Aciv 



A 



1. 



Mgll 



(4) 



The other 'hybrid' approach, which maximizes the 
survey path of this Monte Carlo comparison, is to 
introduce a second random quasar (with z < 1.5) 
to cover the redshift path at z < z m i n .GRB in 
the GRB spectrum. In these cases, the mini- 
mum quasar redshift is z min ,Q SO = (1 + z GRB ) x 
^civMmIh ~ 1- Finally, for very high-redshift 
GRBs (zgrb > 2.2) the second quasar has to be 
chosen such that z min ,Qso = min[(l + zgrb) x 
^civ/^Mgn — 1' 2] and z maXi Q S o = (1 + zgrb) x 
•^Lya/^Mgii — 1' which allows us to select at least 
50 QSOs covering the desired redshift path cover- 
age. 

We ran ten thousand Monte Carlo iterations 
using both approaches and we recorded for each 
iteration the number of Mg II absorbers recov- 
ered. We performed this analysis for each of the 
GRB samples. Figure [7]prcscnts our outcomes us- 
ing the hybrid approach, though no relevant dif- 
ferences are present using the truncated redshift 
path. 

The results indicate that the incidence of Mg II 
absorbers detected in our independent Sample I 
are consistent with the results along quasar sight- 
lines. In fact, we recovered a slightly greater num- 
ber of absorbers on average along the quasar sight- 
lines. Furthermore, the analysis shows that there 
is no statistically significant discrepancy between 
the expected total number of absorbers along the 
QSOs and the full parent GRB sample (Sample 
F). In 6% of our simulated quasar lines of sight, we 
observed a number of absorbers equal to or larger 
than the Sample F (corresponding to a 1.6<r sig- 
nificance). Only in the case of the high-resolution 
sub-sample, Sample H is there a statistically sig- 
nificant excess. This sample, however, is domi- 
nated by the sightlines analyzed in previous works 
(e.g. P06). We discuss this result further in the 



following section. A summary of our Monte Carlo 
analysis is given in Table [2] 

It is further illuminating to estimate the statis- 
tical significance of the Mg II enhancement along 
GRB sightlines as a function of historical time. 
Figure [8] shows the results of a Monte Carlo anal- 
ysis for each year, where we include all GRBs from 
that year and any previous. Until the end of 2006 
a significant (> 3tr) excess was present. Since 
that time, the statistical significance has steadily 
declined and the current full sample (which has 
several times the survey path of P06) has only a 
modest statistical significance. At present, we do 
not find a statistically significant difference in the 
incidence of strong Mg II absorbers between GRB 
and quasar sightlines. 

5.3. Other Characteristics of the Mg II 
GRB Sample 

In Fig. [9] we present the cumulative distribu- 
tions of the equivalent widths and relative veloc- 
ities for the strong Mg II full sample. The lat- 
ter is calculated assuming that these interven- 
ing systems are local to the GRB environments 
and are moving at such velocity towards the ob- 
server to mimic a low er redshift system (see also 
Cucchiara et al. 20091 for the intrinsic properties 
of a small sample of such systems). As pre- 
viously observed, more than 50% of the inter- 
vening systems would require ejection velocities 
larger than 50,000 km s _1 , making very unlikely 
an intrinsic origin of these absorbers. Recently, 
Bergeron et al.l (|2011l ). based on similar distribu- 
tion of strong Mg II absorbers along blazars, have 
suggested a possible theoretical model for produc- 
ing such high relative velocities. 

The red curves represent similar quantities from 
our Monte Carlo analysis of the QSO sightlines. A 
KS-test analysis shows for both metrics that the 
GRB and QSO absorbers are consistent with hav- 
ing been drawn from the same parent population 
(P KS = 0.48 and P KS = 0.39, for the W r and the 
projected velocity respectively). 

6. Discussion and Conclusion 

We have presented the largest compilation to 
date of GRB spectroscopic data, more than one 
hundred spectra including data from previous 
published works, proprietary datasets, and pub- 
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licly available datasets not yet published. We 
have leveraged this dataset to investigate the puz- 
zling excess of strong Mg H absorbers along GRB 
sightlines as first noted bv lProchter et al.l(|2006bh . 
Most importantly, we have performed such analy- 
sis on a fully independent dataset to the original 
P06 study in order to test their findings. 

This independent sample, our Sample I, com- 
prises 83 GRB lines of sights, yielding a red- 
shift path length Az — 44.9 over the interval 
z = 0.4 — 2.2. Along these spectra, we detect 
only 8 absorbers, for a total incidence of strong 
Mgll absorbers (W r > lA) oit{z),i = 0.18. This 
incidence lies in good agreement with estimations 
alo ng QSO lines of s ight ta ken from the latest work 
bv IZhu fc Menardl (120121) {l{z) QSO = 0.26). No 
excess has been identified in the independent sam- 
ple and, therefore, we do not confirm the original 
findings of P06 that an excess of Mg II absorbers 
lie along GRB sightlines. 

It is likely that the earlier works on the inci- 
dence of Mg II absorption along GRB sightlines 
were biased by a remarkable, statistical fluke. In 
particular, the presence of a small set of lines 
of sight with multiple absorbers app ears to have 
drive n the results (as suggested by iKann et aT 
2010h . 



Even including the original P06 data (i.e. our 
full dataset, Sample F), which maximizes the red- 
shift path coverage observed along GRBs sight- 
lines (Az = 55.5), we estimate i(z) : p = 0.36 ± 
0.09, which corresponds to an excess of strong 
Mg II absorbers by a factor ~ 1.5 over QSO sight- 
lines (at 90% c.l. for Poisson distribution). We 
tested the significance of this excess using a Monte 
Carlo analysis and find that 6% of random QSO 
samples exhibit as many absorbers as the GRB 
survey. This suggests the null hypothesis is ruled 
out at < 2cr confidence level. In conclusion, the 
data no longer demand a different incidence of 
strong Mg II absorption along GRB and QSO 
sightlines. 

We wish to emphasize that the P06 analysis was 
not inherently flawed. Indeed, if we restrict our 
analysis to the set of high-resolution data, which 
has large overlap with the P06 sample, we find a 
significant excess (« 3 times) at a high statistical 
significance (« 4<r). At face value, this could sug- 
gest that we have underestimated l{z) for the low- 
resolution sample, e.g. because we mis-estimated 



our sensitivity to lA absorbers. Our sample of 
low-resolution data, however, includes a large di- 
versity of S/N. In order to investigate the effect of 
these diversity on our detection rate, we degraded 
the spectra in Sample H to the lowest S/N and res- 
olution for which we are able to estimate the red- 
shift path length (e.g. the ALFOSC spectrum of 
GRB 050802, which has S/N=7 and R « 440): all 
the strong Mg n doublets could still be detected 
at 5ct level. Furthermore, we have identified many 
additional Mg II absorbers in these spectra (Ta- 
ble [5]) where the selection criteria are not fully sat- 
isfied. We also established our completeness level 
and the reliability of our automatic searching al- 
gorithm creating a larger (^5000) set of spectra, 
derived by the original full sample, where we ran- 
domly injected mocked doublet profiles with dif- 
ferent equivalent widths. The automatic identi- 
fication process recovered ~ 98% of the mocked 
features. At this stage, we suspect that the few 
lines of sight observed with high-resolution spec- 
trographs were simply "peculiar" with respect the 
presence of strong Mg II doublets. Surely a larger 
collection of such data (e.g. the sample building 
with X-Shooter) will allow for an independent test 
of the high-resolution results. 

It is also worth noting, in this context, that 
other authors have explored whether the bright- 
ness of the GRB afterglow correlates with the 
presence of intervening Mg II absorption, i.e. 
to bias the observ ations towards such sightlines. 
Kann et al.l (|201dh have investigated the optical 



properties of these GRBs in relation to the pres- 
ence/absence of Mg II absorbers and the pos- 
sibility that GRB optical afterglows brightness 
may be boosted due to gravitational len sing (see 
IPorciani et aIll2007bl:lMenard et al.ll2008l ). In par- 
ticular they compared the absolute mean i?-band 
magnitude (estimated at one day post-burst and 
normalized at z = 1) of GRB with strong ab- 
sorbers and without (which usually present weak 
absorbers). For this purpose they used afterglow 
spectra obtained with echelle spectrographs which 
provide high-S/N. No appreciable difference was 
noticed betw een the two sam ples. While we defer 
the reader to M enar dl (|2005h for a quantitative es- 
timate of possible gravitational lensing effects, we 
note that considering only the lines of sight with 
strong abso rbers, our S a mple H extends the origi- 
nal work of Kann et al.l ( 2010l ) by only one object, 
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leading to inconclusive progress on this aspect due 
small size samples. 

Moreover, we compared the equivalent width 
distribution of the detected absorbers in our 
Sample F and our Monte Carlo analysis: a 
Kolmogorov-Smirnov test shows that no signifi- 
cant difference is present between the two samples 
(Pks — 0.48). Similarly, if considering the rela- 
tive velocity of the two populations of absorbers 
as they were, instead of intervening, moving at 
high velocity towards the observer so as to mimic 
a lower redshift we also do not find any particular 
difference (Pks = 0.39), further disfavouring an 
intrinsic nature for the absorbers. 

Undoubtedly, the most robust results arc ob- 
tained from high S/N, high resolution (Echelle or 
Echellette) data, of which we only have a limited 
sample for GRB afterglows to-date. For this rea- 
son new samples (such as that being gathered by 
X-shooter) obtained at high resolution will provide 
an important test of our conclusions. 
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B.Zhu for the useful comments and to have pro- 
vided the best to date compilation of high signal 
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Table 1 

List of objects considered for the Mg n analysis 



GRB 


Z GRB 


Telescope 


Instrument 


Resolution 

(A) 


S/N a 


Rclerence 


111229A 


1.380 


Gemini 


CMOS 


5.8 


6.7 


this work 


111107A 


2.893 


Gemini 


CMOS 


5.8 


3.5 


this work 


111008A 


4.989 


Gemini 


CMOS 


5.8 


3.8 


this work 


110918A 


0.982 


Gemini 


GMOS 


5.8 


20 


this work 


110731A 


2.830 


Gemini 


CMOS 


5.8 


26 


this work 


110726A 


1.036 


Gemini 


GMOS 


5.8 


9 


this work 


110213B 


1.083 


Gemini 


GMOS 


3.4 


5 


this work 


110213A 


1.460 


Bok 


FAST 


6 


20 


(4) 


110205A 


2.214 


Lick 


KAST 


11 


14 


(4) 


100906A 


1.727 


Gemini 


GMOS 


5.8 


21 


this work 


100901A 


1.408 


Gemini 


GMOS 


3.4 


6 


this work 


100814A 


1.438 


MAGELLAN 


MagE 


1.8 


10 


this work 


100513A 


4.798 


Gemini 


GMOS 


5.8 


17 


this work 


100418A 


0.624 


VLT 


X-Shooter 


0.86/0.72/2* 


12 - 38 


(19) 


100414A 


1.368 


Gemini 


GMOS 


5.8 


11 


this work 


100302A 


4.813 


Gemini 


GMOS 


5.8 


3 


this work 


100219A 


4.667 


Gemini 


GMOS 


1.6 


1.2 


this work 


091208B 


1.063 


Gemini 


GMOS 


5.8 


26 


(2) 


091109A 


3.076 


VLT 


FORS2 


13 


3 


this work 


091029 


2.752 


Gemini 


GMOS 


5.8 


39 


this work 


091024 


1.092 


Gemini 


GMOS 


5.8 


55 


(2) 


091020A 


1.713 


NOT 


ALFOSC 


13 


7 


this work 


090926A 


2.106 


VLT 


X-Shooter 


1.0 


15 - 30 


(6) 


090902B 


1.822 


Gemini 


GMOS 


4 


14 


this work 


090812A 


2.454 


VLT 


FORS2 


13 


15 


(18) 


090529A 


2.625 


VLT 


FORS2 


13 


4 


(18) 


090519A 


3.851 


VLT 


FORS2 


13 


3 


(18) 


090516A 


4.109 


VLT 


FORS2 


13 


24 


(18) 


090426 


2.609 


Keck 


LRIS 


5.5 


8 


this work 


090424 


0.544 


Gemini 


GMOS 


5.8 


22 


this work 


090323 


3.567 


Gemini 


GMOS 


5.8 


17 


this work 


090313 


3.375 


Gemini 


GMOS 


5.8 


11 


this work 


081222 


2.771 


Gemini 


GMOS 


5.8 


21 


this work 


081029 


3.847 


Gemini 


GMOS 


5.8 


42 


(2) 


081008 


1.967 


Gemini 


GMOS 


3.4 


35 


(2) 


081007 


0.529 


Gemini 


GMOS 


5.8 


31 


(2) 


080928 


1.690 


Gemini/VLT 


GMOS/FORS2 


5.8/13 


8/25 


(2)/(3) 


080916A 


0.689 


VLT 


FORS1 


13 


5 


(18) 


080913A 


6.700 


VLT 


FORS2 


13 


2.5 


(13) 


080905B 


2.374 


VLT 


FORS1 


13 


13 


(18) 


080810 


3.350 


Keck 


HIRES 


0.18 


16 


(18) 


080805 


1.505 


VLT 


FORS2 


13 


3 


(3) 


080804 


2.205 


Gemini 


GMOS 


5.8 


17 


(2) 


080721 


2.608 


TNG 


Dolores 


8.1 


9 


(18) 


080710 


0.845 


Gemini 


GMOS 


3.4 


35 


(2) 


080707 


1.234 


VLT 


FORS1 


13 


5 


(3) 


080607 


3.036 


Keck 


LRIS 


4 


11 


(3) 


080605A 


1.639 


VLT 


FORS2 


13 


30 


(3) 


080604 


1.416 


Gemini 


GMOS 


5.8 


4 


(2) 


080603B 


2.686 


NOT 


ALFOSC 


13 


41 


(3) 
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Table 1 — Continued 



GRB 


ZGRB 


Telescope 


Instrument 


Resolution 

(A) 


S/N a 


Reference 


080603A 


1.688 


Gemini 


GMOS 


5.8 


38 


(2) 


080520 


1.545 


VLT 


FORS2 


13 


5 


(3) 


080413B 


1.100 


Gemini 


GMOS 


3.4 


2 


(2) 


080413A 


2.433 


Gemini 


GMOS 


4 


14 


(2)/(9) 


080411 


1.030 


VLT 


FORS1 


13 


60 


(18) 


080330 


1.513 


NOT 


ALFOSC 


13 


18 


(3) 


080319C 


1.949 


Gemini 


GMOS 


3.4 


4 


(2) 


080319B 


0.937 


Gemini/ VLT 


GMOS/UVES 


5.8/0.13 


45/70 


(2)/(9) 


080310A 


2.4272 


VLT 


UVES 


0.13 


15 


(9) 


080210 


2.6419 


VLT 


FORS2 


13 


33 


(3) 


071122 


1.141 


Gemini 


GMOS 


5.8 


12 


(2) 


071117 


1.334 


VLT 


FORS1 


13 


4 


(3) 


071112C 


0.823 


Gemini 


GMOS 


4 


3 


(2) 


071031 


2.692 


VLT 


UVES/FORS2 


0.13/13 


70/40 


(3) 


071020 


2.145 


VLT 


FORS2 


13 


6 


(3) 


071010B 


0.947 


Gemini 


GMOS 


5.8 


15 


(2) 


071003 


1.6044 


Keck 


LRIS 


5 


34 


(20) 


070810A 


2.170 


Keck 


LRIS 


5 


6 


(21) 


070802 


2.453 


VLT 


FORS2 


13 


8 


(3) 


070721B 


3.626 


VLT 


FORS2 


13 


6 


(18) 


070611 


2.039 


VLT 


FORS2 


13 


15 


(3) 


070529 


2.498 


Gemini 


GMOS 


3.4 


12 


(2) 


070506 


2.306 


VLT 


FORS1 


13 


3 


(3) 


070411 


2.954 


VLT 


FORS1 


13 


6.5 


(3) 


070318 


0.836 


Gemini 


GMOS 


4 


8 


(2) 


070306 


1.496 


VLT 


FORS2 


13 


4 


(3) 


070125 


1.547 


Gemini/ VLT 


GMOS/FORS1 


5.8/ 


6/15 


(10)/(3) 


070110 


2.352 


VLT 


FORS2 


13 


30 


(3) 


061121 


1.314 


Keck 


LRIS 


13 


28 


(3) 


061110B 


3.434 


VLT 


FORS1 


13 


11 


(3) 


061110A 


0.758 


VLT 


FORS1 


13 


6 


(3) 


061007 


1.261 


VLT 


FORS1 


13 


6 


(3) 


060927 


5.468 


VLT 


FORS1 


8 


2.5 


(1) 


060926 


3.205 


VLT 


FORS1 


13 


8 


(1) 


060908 


1.884 


Gemini 


GMOS 


5.8 


8 


this work 


060906 


3.686 


VLT 


FORS1 


13 


5.8 


(1) 


060904B 


0.703 


VLT 


FORS1 


13 


12 


(18) 


060729 


0.543 


Gemini 


GMOS 


3.4 


26 


(2) 


060714 


2.711 


VLT 


FORS1 


13 


50 


(3) 


060708 


1.923 


VLT 


FORS2 


13 


5 


(3) 


060707 


3.425 


VLT 


FORS2 


13 


7 


(3) 


060607A 


3.047 


VLT 


UVES 


0.13 


43 


(9) 


060526 


3.221 


VLT 


FORS1 


13 


38 


(1) 


060522 


5.111 


Keck 


LRIS 


5 


2.3 


(1) 


060512 


2.092 


VLT 


FORS1 


13 


3 


(3) 


060510B 


4.922 


Gemini 


GMOS 


5.8 


8.4 


(2) 


060502A 


1.515 


Gemini 


GMOS 


5.8 


8 


(2) 


060418 


1.489 


Gemini/ VLT 


GMOS/UVES 


5.8/0.13 


86/60 


(2)/(H) 


060210 


3.912 


Gemini 


GMOS 


5.7 


26 


(2) 


060206 


4.046 


NOT 


ALFOSC 


13 


40 


(3) 
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Table 1 — Continued 



GRB 


Z GRB 


Telescope 


Instrument 


Resolution 

(A) 


S/N a 


Reference 


060124 


2.296 


Keck 


ESI 


13 


8 


(3) 


060115 


3.5328 


VLT 


FORS1 


13 


10 


(3) 


051111 


1.5489 


Keck 


HIRES 


0.18 


20 


(12) 


050922C 


2.1996 


VLT 


UVES 


0.13 


12 


(9) 


uuuyuo 


o.ooy 


Gemini /Keck 






Q /I 9 


( 2) / this work 


050820 


2.614 


VLT 


UVES 


0.13 


23 


(9) 


050802 


1.711 


NOT 


ALFOSC 


13 


7 


(3) 


050801 


1.559 


Keck 


LRIS 


5 


5 


(3) 


050730 


3.9687 


VLT 


UVES 


0.13 


40 


(9) 


050401 


2.896 


VLT 


FORS2 


13 


23 


(3) 


050319 


3.240 


NOT 


ALFOSC 


13 


6 


(3) 


030429 


2.655 


VLT 


FORS1 


13 


7 


(1) 


030323 


3.372 


VLT 


FORS1 


13 


8 


(13) 


030226 


1.986 


VLT 


FORS1 


13 


30 


(14) 


021004 


2.323 


VLT 


UVES 


0.13 


10 


(9) 


020813 


1.255 


VLT 


UVES 


0.13 


60 


(15) 


010222 


1.477 


Keck 


ESI 


0.6 


4 


(16) 


000926 


2.038 


Keck 


ESI 


0.6 


12 


(17) 



a Signal-to-noisc ratio is estimated as the median at the continuum level over a wavelength range 
clean of telluric lines 



References. — (1) Jakobsson ct al N2006T ); f2) ICucchiaralll2010l); f3HFvnbo et al.H2009l) and ref- 
erence ther ein; (4) ICucchiara et al.l feoilaf); (5)lde Ugarte Postigo et al l ll201lt); (6 ) D'Elia et all 
j20ld): (7lD'EUa, et alJ ll201lh: (8) iThoene et al l ll200Sfl; (91 IVergani et i al] 1120091): (TO) 
ICenko et al.H200g|); (TDlVreeswiik et alj<2007l); (T2)|Prochaska et al.| J2007T): (T3HVreeswiik et al.l 



(14) 



Klosc ct al. 



<2004ft; 
ICastro 

f20) |Perle" et al 



ICastro et ail J2003T); (18 ) lde Ugarte Postigo et al.l 
I j2008T) ;f21') lMilvang-Jensen et al 



(15) iBarth et alj ll2003ft ; ( 16) iMirabal et al.l J2002 
' ' 2012a|) ;7l9) Ide Ugarte Postigo et al.1 



(2012) 



Table 2 
Studied Sample 







Number a 
of GRBs 


Az ■ 

l.OA 


N ■ 

l.OA 






Sample 


1 


83 


44.9 


8 


0.18 ±0.06 


0.26 


Sample 


F 


95 


55.5 


20 


0.36 ± 0.09 


0.24 


Sample 


11 


18 


20.3 


13 


0.64 ±0.26 


0.25 


Sample 


L 


79 


35.3 


7 


0.19 ±0.08 


0.12 



a The total number of GRB lines of sight in each sample corresponds 
to the sum of all those GRBs where £(z) 7^ 0. 



Note. — Summary of our Mg n search: the sample name and the num- 
ber of lines of sight included are listed in the first two columns; the total 
redshift path density explored and the number of absorbers identified are 
listed in the third and forth column. Based on these we could determine 
the incidence of the absorbers in each sample and compare it with the 
expected incidence along our QSOs sample (last two columns). 



17 



Table 3 

List of excluded lines and regions in the redshift pathlength estimate 



Description 


Arest (A) 


Description 


Arest(A) 


Nv 


1238,1242 


C I 


1560, 


Sn 


1250,1253,1259 


Fc II 


1608,1611,2249,2260,2344,2374,2382,2586,2600 


Si II 


1260,1304,1526,1808 


Al ii 


1670 


Si ii* 


1264,1309,1533,1816 


Al in 


1854,862 


Oi 


1302 


Cr ii 


2017,2026,2056,2066 


Ni ii 


1317,1370,1454,1703,1709,1741,1751 


Zn II 


2026,2062 


C ii 


1334 


Ni II* 


2217 


C ii* 


1335 


Mnll 


2576,2594,2606 


Si iv 


1393,1402 


Band B a 


6860 - 7000 


C iv 


1548,1550 


Band A a 


7600 - 7704 


Atm. Band" 


8130 - 8323 


Atm. Band" 


8930 - 9020 



a Atmospheric absorption bands from the HIRES telluric line list: http: //www2 . keck.hawaii . edu/insthires/makeewww/Atmosphere/atmabs . txt. 
The indicated wavelengths arc (obviously) independent of the GRB redshift. 



Table 4 

List of Redshifts Intervals where g(z) = 1 for lA Mg n survey 



GRB z G rb 


Zstart 


z end 


030226 1.986 


1.30132 


1.31706 




1.32636 


1.32922 




1.35426 


1.38430 




1.41005 


1.41935 




1.42936 


1.45655 




1.46656 


1.48158 




1.50161 


1.54667 




1.55669 


1.62393 




1.64038 


1.64682 




1.66971 


1.71119 




1.72407 


1.76127 




1.77200 


1.77700 




1.78773 


1.81134 




1.82922 


1.85211 




1.86284 


1.86356 




1.87428 


1.92364 




1.94367 


1.97300 




1.99231 


2.14664 




2.16738 


2.18741 




2.20959 


2.35908 




2.37124 


2.39413 




2.41774 


2.44850 




2.48927 


2.48999 




2.50644 


2.52718 




2.54793 


2.70959 




2.75894 


2.75966 




2.78040 


2.78112 




2.78684 


2.90057 
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Table 5 
Intervening systems 



GRB 


Z GRB 


Zabs 


iy r (2796) 


W r (2803) 


Statistical 


Other 








(A) 


(A) 


sample 


transition 


010222 


1.477 


1.156 


2.22(0.14) 


1.69(0.11) 


F,H 


Fe II 


020813 


1.255 


1.224 


1.58(0.03) 


1.43(0.03) 


F,H 


Mg i,Fe II 


021004 c 


2.3295 


0.555 


0.66(0.045) 


0.36(0.034) 


N 


Mg l,Fc II 






1.380 


1.637(0.020) 


1.574(0.043) 


F,H 


Mg i,Mnll,Fe n 






1.6026 


1.407(0.024) 


1.02(0.013) 


F,H 


Mg I, Fe ll,Mnll 


030226 d 


1.986 


1.043 


0.68(0.25) 


0.41(0.25) 


N 


Al n 






1.963 


2.22(0.10) 


2.47(0.10) 


N 


Mg i, C iv, Si II 


050730 c 


3.9687 


1.7732 


0.927(0.030) 


0.718(0.016) 


N 


Mg l,Fe II 






2.2531 


0.783(0.650)° 


0.677(0.017) 


N 


Si n, Al n, Fe n, Mg I 


050820 c 


2.6147 


0.6915 


2.723(0.007) 


1.576(0.031) 


F,H 


Mg I 






1.4288 


1.203(0.023) 


1.265(0.026) 


F,H 


Mg i,Fe II, Al in 






1.6204 


0.277(0.024) 


0.214(0.008) 


N 


Mg i,Fe II, Zn II, Si n 






2.3598 


0.424(0.306)° 


0.517(0.024) 


N 


Fe II, Si II, Zn n,C iv 


050908 


3.339 


1.548 


1.21(0.02) 


0.92(0.02) 


F,H 


Fe II 


050922C C 


2.1996 


0.6369 


0.187(0.018) 


0.121(0.011) 


N 


Mg i,Fe n 






1.1076 


0.476(0.029) 


0.422(0.19) 


N 


Mg i,Fe II 






1.5670 


0.121(0.080)° 


0.088(0.007) 


N 


C iv,Fe II 


051111 


1.55 


1.190 


1.56(0.02) 


1.92(0.01) 


F,H 


Mg i,Fe II 






0.827 


0.39(0.02) 


0.29(0.01) 


N 


Mg I 


060418'' 


1.489 


1.107 


1.84(0.2) 


1.58(0.1) 


F,H 


Mg i,Fe II, Zn n, Al m,Al n 






0.6559 


1.52(0.3) 


2.15(0.4) 


F,H 


Fe II 






0.603 


1.49(0.2) 


1.47(0.1) 


F,H 


Fe II 


060502A 


1.515 


1.147 


2.39(0.12) 


2.87(0.12) 


F,L,I 


Mg I 






1.078 


0.61(0.12) 


0.49(0.12) 


N 








1.044 


1.90(0.15) 


1.92(0.16) 


F,L,I 


Fe l,Mnll,Mg I 


060607A C 


3.0748 


1.5103 


0.124(0.011) 


0.144(0.007) 


N 


Fe II 






1.8033 


1.916(0.006) 


1.600(0.015) 


F,H,I 


Mg i,Fc n, Al in 






2.2783 


0.210(0.058) 


0.298(0.013) 


N 


Fe n,Al m,Al II, C iv, Si n,Si iv 


060906 


3.685 


1.2659 


1.63(0.28)° 


1.63(0.28)° 


N 


Mg I 


060926 


3.2 


0.924 


2.49(0.62)° 


2.49(0.62)° 


N 


Mg i,Fc I 






1.7954 


3.27(0.69) 


3.71(0.87) 


N 


Mg l,Fe n,Mnll 






1.8289 


1.27(0.11) 


0.72(0.07) 


N 


Mg I 


061007 


1.261 


1.065 


3.14(0.53)° 


4.48(0.65) f ' 


N 


Mg l,Fe n,Mnll 


070529 


2.498 


1.414 


0.20(0.02) 


0.09(0.02) 


N 




070506 


2.306 


1.600 


1.92(0.04) 


1.65(0.05) 


N 


Al in 


070611 


2.039 


1.297 


2.65(0.27) 


1.99(0.23) 


N 


Mg i,Fe n 


070802 


2.45 


2.0785 


0.82(0.12) 


0.82(0.12) 


N 


Al n, Ni ii, Mg i,Fe n 






2.2921 


0.55(0.15) 


0.55(0.22) 


N 


Ni n, Al ni,Cr n, Fc n 


071003 


1.604 


0.372 


2.28(0.19) 


1.91(0.19) 


F,L,I 


Mg I 






0.943 


0.61(0.05)° 


0.36(0.05) 


N 


Mg I 






1.101 


0.80(0.06) 


0.64(0.05) 


N 


Mg I 


071031 c 


2.6922 


1.0743 


0.330(0.016) 


0.206 (0.008) 


N 


Fe ii 






1.6419 


0.806(0.014) 


0.586 (0.052) 


N 


Fe ii, Al iii,C iv 






1.9520 


0.743(0.016) 


0.612 (0.016) 


N 


Mg i,Fe ii 


080310 c 


2.4272 


1.6711 


0.421(0.012) 


0.366 (0.016) 


N 


Mg i,Fe ii, Al ii, Si II, C iv 


080319B C 


0.9378 


0.5308 


0.614(0.001) 


0.350 (0.002) 


N 


Mg i,Fc ii 






0.5662 


0.083(0.003) 


0.029 (0.001) 


N 


Mg i,Fe ii 






0.7154 


1.482(0.001) 


0.736 (0.003) 


F,H,I 


Mg i,Fe ii 






0.7608 


0.108(0.002) 


0.039 (0.002) 


N 


Fe II 


080319C 


1.95 


0.8104 


2.04(0.52) 


1.64(0.42) 


N 


Fe ii,Mno 
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Table 5 — Continued 



GRB 


ZGRB 


^abs 


W r (2796) 


Wr(2803) 


Statistical 


Other 








(A) 


(A) 


sample 


transition 


080603A 


1.688 


1.271 


3.11(0.11) 


3.17(0. 13) b 


F,L,I 


Mg l,Fe II 






1.563 


0.77(0.01) 


0.92(0.01) 


N 


Fe I 


080605 


1.64 


1.2987 


1.08(0.11) 


0.77(0.10) 


F,L,I 


Fc II 


080607A 


3.036 


1.341 


3.0(0.08) 


1.26(0.05) 


F,L,I 


Mg I 


080805A 


1.505 


1.197 


8.2(0.92)° 


8.2(0.92)" 


N 


Mnii,Fe n 


080905B 


2.374 


0.618 


6.65(0.2)" 


6.65(0.1)" 


N 


Mg I 


080928 


1.691 


0.736 


9.54(0.25)" 


9.54(0.25)" 


N 


Mg i,Fe II 


081222 


2.77 


0.8168 


0.52(0.01) 


0.28(0.11) 


N 


Mg i,Fe II 






1.0708 


1.46(0.23) 


0.61(0.21) 


F,L,I 


Fe II 


091208B 


1.063 


0.784 


0.65(0.43) 


1.03(0.43) 


N 


Mg I 


100814A 


1.44 


1.1574 


0.426(0.04) 


0.379(0.04) 


N 


Mg I 


100901A 


1.408 


1.314 


1.74(0.17) 6 


i.53(o.i6) b 


N 


Fe n, Mg I 


100906A 


1.64 


0.994 


0.87(0.1) 


l.ig^o.i) 6 


N 




110918A 


0.982 


0.877 


2.65(0.20) 


2.82(0.20) 


N 


Mg i,Fe II 



a Equivalent Width measurement is largely effected by blending. For these lines we 
report the total EW for the doublet. 

b Equivalent Width measurement is lightly effected by blending. EW values are derived 
via deblending procedure using gaussian fit of the two lines (either the other member of 
the doublet or other lines) via the IRAF splot tool. 

C UVES 

d Also see lKlose et""aH i2004bT) 
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Fig. 1. — Comparison between 4 different GRB spectra obtained with different spectrographs and different 
resolving power. From top to bottom: GRB 050922C observed with ALFOSC; GRB 060210 observed with 
Gemini/GMOS; GRB 100418A observed with the UV arm of VLT/X-Shooter; section of GRB 060607A 
observed with VLT/UVES. In all the panes the grey curve represent the associated la error spectrum. 
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GRB 060210 (z = 3.912) 




6000 



6500 7000 

Observed wavelength 



7500 



8000 



Fig. 2. — GRB 060210 Gemini spectrum. This example shows our excluded regions for the purpose of 
estimating the survey path for intervening Mg II absor ption taking into account the host absorption lines 
(in gold, from the tabulation of IChristensen et al. I (Son!)) as well as telluric lines (in gray) . The red curve is 
the la spectrum associated with the data. 
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Fig. 3. — Le/t:Rcdshift path density for 1 A rest-frame equivalent widths at 5cr detection limit for the Sample 
F (blue) and Sample I (black). Dotted vertical lines represent the quasar selection regions, where the Mg II 
doublet is detectable in the SDSS spectral coverage. Right: similar plot for Sample H and Sample L. 
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Recovered doublets 
Unrecovered doublets 



V300 AlfoscDolores Kast R400 B600 V1200 MagE ESI HiRes UVES 




V300 AlfoscDolores Kast 



B600 V1200 MagE 



HiRes UVES 



Fig. 4. — Top: Number of recovered (black) and unrecovered (red) strong doublets (W r > lA) in our 
mock sample ordered by instrument resolution. Most of the missed doublets, besides being in searchable 
regions of the spectra (g(z) — 1), are missed or misidentified due to self-blending with other features (like 
other intervening systems metal lines, or wings of GRB host features). Therefore the automatic procedure 
usually fails to identify both members of the doublet due to the low-resolution of the instrument. Bottom: 
Completeness level ordered by spectral resolution. From our mock sample we derived a final completeness 
level of - 98%. 
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Fig. 5. — £(z) evolution of intervening Mg II absorbers (W2796 > 1A) for our sample of GRB sightlines: 
triangles and square symbols refer to the Sample I and Sample F respectively. The red curve show s the 
evolution of the Mg 11 incidence along quasar sightlines as recently computed bv lZhu fc Menardl (|2012t) . We 
derive an average £(z) = 0.20 for Sample I, in agreement with the prediction, while t{z) — 0.36 for Sample 
F, indicating a slight overabundance of absorbers compared to the QSOs. 
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J 1 i - J j. | | | | | | | | | | | | | | | L 



0.0 0.5 1.0 1.5 2.0 

Redshift 

Fig. 6. — Cumulative distribution of strong Mg n absorbers along GRB sightlines for Sample I and Sample F 
(black and blue solid curves, respectively). These are compared to the predicted incidence based on measure- 
ment along QSO lines of sight (dashed curves). The independent Sample I actually shows fewer absorbers 
than expected while a modest excess remains in Sample F. Neither result corresponds to a statistically 
significant difference from the QSO results. 



2G 



1200 
| 1000 

g 800 

• rH 

5 600 

3 400 

| 200 


g 1200 
1 1000 
800 

o 

1 600 

| 400 

| 200 





i 



J 



H 



10 20 30 
Number of absorbers 



5 10 15 20 25 30 
Number of absorbers 



Fig. 7. — The black curves show the distribution of recovered Mg n absorbers along Monte Carlo realizations 
of quasar sightlines designed to match the g(z) survey path of the GRB samples (from top left clockwise, 
Sample F, Sample I, Sample L, and Sample H). The red-dashed lines trace the mean number of absorbers 
for each distribution and the shaded regions represent the 1,2, and 3a confidence interval assuming Poisson 
statistics. The solid arrow in each panel denotes the number of Mg II absorbers detected for each subsample 
of GRBs. Only the high-resolution Sample H exhibits a statistically significant excess, but we caution that 
this sample has substantial overlap with the original P06 work. 
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Included Year 



2010 2012 



Fig. 8. — Top : Observed excess in the incidence of strong Mg II absorption along GRB sightlines relative 
to that predicted from observations along quasar sightlines. This is shown as a function of historical time 
where each bin includes all the GRB lines of sights until December 31 s * of the specified year, as extracted 
from Sample F. The filled star marks the results published by P06. Bottom : Confidence level at which the 
excess factor has been detected based on the Monte Carlo analysis described in § 15.21 
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Fig. 9. — Left : Rest-frame cumulative equivalent width distribution of the Mg n absorbers in Sample F 
(black) and the QSOs absorbers from our Monte Carlo analysis (red). Right : Cumulative distribution of 
the relative velocity (black for the GRB and red for the QSO absorbers), assuming every absorber is local to 
the QSO or the GRB host galaxy and is moving towards the observer mimicking a foreground intervening 
system at lower redshift). Kolmogorov-Smirnov tests reveal that the properties of the Mg II absorbers along 
GRB and QSO sightlines are consistent with having been drawn from the same parent population. 
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Fig. 10a. — Line profiles of strong Mg n transitions along GRB sightlincs. 
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Fig. 10b.- 
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Fig. 9c- 
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Fig. lOd. 
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Fig. lOe. 
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Fig. 10f.- 
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Fig. lOg, 
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